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Abstract 

We present a study of the branching fraction of the decay B° — > p°p° and the fraction of 
longitudinally polarized p° mesons in this decay. The results are obtained from the final data 
sample containing 772 x 10 6 BB pairs collected at the Y(4S) resonance with the Belle detector at 
the KEKB asymmetric-energy e + e~ collider. We find 166 ±59 B° — > p°p° events, corresponding to 
a branching fraction of B(B° ->• p°p°) = (1.02±0.30 (stat)±0.22 (syst))xl0~ 6 with a longitudinally 
polarization fraction f^ = 0.21^Q22 (stat) ± 0.11 (syst). The branching fraction's upper limit is 

B{B° -> p°p°) < 1.5 x 1(T 6 

at 90% confidence level. We use the longitudinal polarization fraction to determine the Cabibbo- 
Kobayashi-Maskawa matrix angle 4>2 = (91-0 db 7.2)° through an isospin analysis in the B — )• pp 
system. We furthermore obtain 

B(B° -> /op ) x B(f -> 7T+7T-) = (0.86 ± 0.27 (stat) ± 0.15 (syst)) x 10~ 6 , 

corresponding to 149 ± 49 B° — > fop events with a significance of 3.0 standard deviations and 
upper limits at 90% confidence level on the (product) branching fractions for possible decays with 
the same final state, 

B(B° -> vr+Tr^vr+vr") < 11.7 x 10~ 6 , 
B(B° -»• p°7r + 7r-) < 12.2 x 10~ 6 , 
B(B° -»• /ott+tt-) x B(f -> tt + tt~) < 3.1 x 10~ 6 and 
B(B° -> / /o) x B(f -)• tt+tt-) 2 < 0.2 x HT 6 . 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 
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I. INTRODUCTION 



CP violation in the standard model (SM) is due to an irreducible complex phase in 
the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix [l], 01. Mixing-induced CP 
violation in the B sector has been clearly observed by the Belle [!, H[ and BaBar jE], |6[ 
collaborations in the b — > ccs induced decay B° — > J/tpKg, while many other modes provide 
additional information on CP violating parameters. 

Decays that proceed predominantly through the b — > uud transition are sensitive to one 
of the angles of the unitarity triangle, 02 (or a) = arg[(— V t dV t * b ) / (y u dV* b )}. The Belle and 
BaBar collaborations have reported time-dependent CP asymmetries in these modes that 
include decays such as B° vr+vr" p ± ?r T P + P~ 0, EH and af^ ^M±M- A 

feature common to these measurements is that possible loop contributions, in addition to the 
leading order tree amplitude, can shift the measured angle to 02 ff = 02 + A0 2 . Fortunately, 
this inconvenience can be overcome with bounds on A0 2 determined using either an isospin 
analysis [l6| or SU(3) flavor symmetry 17 . 

This analysis is concerned with the branching fraction of B° — > p°p° decays, shown in 
Fig. [H the fraction of longitudinal polarization in these decays and decays into four charged 
pion final states as the p° decays dominantly into two charged pions. 

Since the dominant tree process is color-suppressed, B° — > p°p° is expected to be rarer 
than its isospin partners, making the isospin analysis less ambiguous. The vector state p°p° 
is not a pure CP eigenstate, but rather a superposition of CP-even and -odd states, or three 
helicity amplitudes, with only the longitudinal one being a pure CP eigenstate. In general, 
the different helicity amplitudes can be separated through an angular analysis. Here, we 
use the helicity basis where the angle Hk \ k=lt2 , defined as the angle between the 7r + 's and 
the .B's flight directions in the rest frame of the k th p° (Fig. [2]), can be used to separate 
longitudinally (CP-even) from transversely (CP-even and -odd) polarized p° mesons. The 
angular distribution is given by 



d 2 T 



Yd cos 6 H id cos 9 H2 



^(1 - f L ) sin 2 9 hi sin 2 9 H 2 + fh cos 2 9 m cos 2 9 



H2 



where fi = |A)| 2 /X] \Ai\ 2 is the fraction of longitudinal polarization; A , A +1 and A_i are 
the helicity amplitudes. 

The SM, using perturbative QCD (pQCD) [l8[ or QCD factorization [ID, 2(|, predicts 
the B° — > p°p° branching fraction to be ~ 1 x 10~ 6 . Analogous to the decay of a B° meson 




FIG. 1: Leading-order tree (left) and penguin (right) diagrams for the decay B° — > p°p°. 
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into two charged p mesons, the decay B° — » p°p° is expected to be predominantly longitu- 
dinally polarized. However, color-suppressed B meson decays into two vector particles are 
especially difficult to predict. Hence, this analysis provides a excellent test of the theoret- 
ical assumptions used in these calculations and improves our understanding of the helicity 
structure of the weak interaction. 

Using a sample containing 657 x 10 6 BB pairs and assuming pure longitudinal polariza- 
tion, the previous Belle analysis set an upper limit on the branching fraction of B° — > p°p° 
decays; B(B° p°p°) < 1.0 x 10~ 6 at 90% confidence level (CL) (2lj._The BaBar collab- 
oration has performed a study of B° — > p°p° decays with 465 x 10 6 BB pairs and found a 
branching fraction B(B° p°p°) = (0.92±0.32(stat) ±0.14(syst)) x 10~ 6 and a longitudinal 
polarization fraction of f L = O^l^^stat) ± 0.04(syst) Q. 

In Sec. m we briefly describe the data set and the Belle detector. The event selection and 
the model used for the branching fraction measurement are described in sections II III and IIV[ 
respectively. Next, the fit result is presented in Sec. [V] followed by validity checks in Sec. IVIl 
where the significance of the fit result is discussed as well. The systematic uncertainties are 
given in Sec. IVIII and a constraint of the CKM phase <p2 is presented in Sec. IVIIfl followed 
by our conclusion. 



II. DATA SET AND BELLE DETECTOR 



This measurement is based on the final data sample containing 772 x 10 6 BB pairs 
collected with the Belle detector at the KEKB asymmetric-energy e + e~ (3.5 on 8 GeV) col- 



lider (23(| . At the T(4S") resonance (a/s = 10.58 GeV), the Lorentz boost of the produced BB 
pairs was /?7 = 0.425 along the z direction, which is opposite the positron beam direction. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprising of CsI(Tl) crystals (ECL) located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K\ mesons and to identify muons 



(KLM). The detector is described in detail elsewhere 24]. 

Two inner detector configurations were used. A 2.0 cm radius beampipe and a 3-layer 
silicon vertex detector (SVD1) was used for the first sample of 152 x 10 & BB pairs, while 
a 1.5 cm radius beampipe, a 4-layer silicon detector (SVD2) and a small-cell inner drift 




FIG. 2: Definition of the helicity angles 6}jk f° r each p, identified by index k = 1,2. 
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TABLE I: Summary of peaking background vetoes. A muon mass hypotheses has been applied to 
specific tracks for the J/ip channel. Here, ^1^2 f° rms the first p° candidate and vr^"7r^ the other. 
X represents any intermediate state or track combination that leads to a 4-body final state. 

Regions vetoed Modes vetoed 

1.85 GeV/c 2 < m(7r+vr-) < 1.89 GeV/c 2 
1.85 GeV/c 2 < m(^+^+) < 1.89 GeV/c 2 
1.85 GeV/c 2 < m^Tr+vr^ ) < 1.89 GeV/c 2 
1.95 GeV/c 2 < m(7r+7r2"7r+) < 1.99 GeV/c 2 
1.95 GeV/c 2 < m(TT^TT+7r^) < 1.99 GeV/c 2 

3.06 GeV/c 2 < m(p + p~) < 3.14 GeV/c 2 
0.478 GeV/c 2 < m(7r 2 7r 3 ) < 0.512 GeV/c 2 



B -> £)+[7r+7r-7r + ]A' 
B -> Bj[7r-7r+7r-]X 

J3 -> i^g[7r+7r-]A 



25 . We use a GEANT- 



chamber were used to record the remaining 620 x 10 6 BB pairs 
based Monte Carlo (MC) simulation to model the response of the detector and determine 
its acceptance [26 ]. 



III. EVENT SELECTION 

We reconstruct B° — > p°p° where p° — > 7r + 7r~. Charged tracks have to fulfill loose 
requirements on the distance of closest approach to the interaction point (IP) of \dz\ < 5.0 cm 
and dr < 0.5 cm along and perpendicular to the z axis, respectively. With information 
obtained from the CDC, ACC and TOF, particle identification (PID) is determined with 
the likelihood ratio Cuj = £.;/ + Cj), where Ci [Cj) is the likelihood that the particle is 
of type % (J). We require £k/it < 0.4, which retains 90% of all pions but only 10% of kaons. 
In addition, we place vetoes on particles consistent with the electron or proton hypotheses. 



SVD requirements of two z hits and one r — hit [27J are imposed on the charged tracks, 
to permit a subsequent measurement of the CP asymmetries. 

Intermediate dipion states are reconstructed above the Kg region with an invariant mass 



0.52 GeV/c 2 < m(7r + 7r~) < 1.15 GeV/c 2 straddling the broad p°(770) resonance [28|. This 
range retains 93% of all p°. Upon combination of two dipion states, a B candidate is formed. 
All remaining particles are associated with the accompanying B meson in the event, referred 
to as -Bxag- 

Reconstructed B candidates are described with two kinematic variables: the beam- 
energy-constrained mass Mb c = y' '(E™^/ c 2 ) 2 — (p^ MS /c) 2 , and the energy difference 
AE = E% MS - -Ebc^m> where E^^ is the beam energy and E% MS Qo§ MS ) is the energy 
(momentum) of the B meson, evaluated in the center-of-mass system (CMS). The B candi- 
dates that satisfy M^ c > 5.27 GeV/c 2 and \ AE\ < 0.1 GeV are selected for further analysis. 

To reduce peaking background coming from charm (b — > c) decays of the B meson with 
a similar final-state topology such as B° — > D~[7t~7t + 7t~]tt + or backgrounds due to particle 
mis-identification, we place vetoes on various combinations of the four charged tracks forming 
our B° candidate, as summarized in Table [B The total efficiency loss due to these vetoes is 
4.4%. The dominant background contribution comes from continuum (e + e~ — > qq, where 
q — u,d, s, c) events. We can use their jet-like topology to separate them from the more 
spherical BB decays using a Fisher discriminant jFs/b constructed from the following seven 
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variables: 



• L2' 11 = J2 cn PcX S cos2 ^p?^ s ,tb; where the sum of the CMS momenta runs over the 
charged tracks (c) and neutral clusters (n) on the tag side; the angle is between the 
particle direction and B thrust direction. 

• cos#tb,to; where the angle is between the B thrust direction and the thrust of the 
tag side. 

• cos 9 b, z , where the angle is between the B flight direction and the z direction. 

• hk° — ^2ij k \Pi\\Pjk\^( cos (^ijk))\^ where pi is the CMS momentum of the i th track from 
the signal side (s), p~j k is the CMS momentum of the j£ h particle from the other side 
(o), 9ij k is the angle between particle i and j k and P 2 is the second order Legendre 
polynomial. For the other side, we distinguish three cases using index k — 0, 1 and 
2, for charged tracks, neutral particles and missing energy (treated as a particle), 
respectively, 

which are shown in Fig. [3J The BB and qq training sample is taken from signal MC 
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FIG. 3: (color online) Simulated (MC) and off-resonance data distributions for the quantities used 
to construct the Fisher discriminant J~s/Bi normalized to have the same area. The solid (blue) 
histograms shows the distribution for BB, while the dashed (red) histograms show distributions 
for (qq) events from data taken below the T(45) resonance. The vertical line implies the cut of 
cos#tb to < 0.9. 



and from data taken below the T(4S) resonance, respectively, and a loose requirement of 
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TABLE II: The reconstruction efficiencies (e) of the signal for the four pion final states calculated 
from Monte Carlo, where the errors are statistical. For B° — > p°p°, we use the fz value obtained 
from the fit to data. 



e (mode) SVD1 (%) SVD2 (%) 

B° p°p° 22.38 ± 0.05 25.38 ± 0.06 

B° -»• /o/o 24.73 ±0.06 28.11 ±0.07 

B° -> /op 23.42 ± 0.05 26.96 ± 0.07 

B° foir+ir- 2.82 ±0.02 3.16 ±0.02 

B° -> p°tt+tt- 2.64 ±0.02 3.02 ±0.02 

B° it + it-it + it- 0.98 ±0.01 1.09 ±0.01 



cos#tb,to < 0.9 is placed that retains 90% signal while rejecting 50% of the continuum 
background. The Fisher discriminant is also required to satisfy —3 < J-s/b < 2, this, 
together with the previous mentioned requirements of \AE\ < 0.1 GeV, M^ c > 5.27 GeV/c 2 
and 0.52 GeV/c 2 < m(7r + 7r~) < 1.15 GeV/c 2 , defines the fit region. 

On average, 1.06 B candidates are reconstructed per signal event. Selecting the best B 
candidate having M^ c nearest the nominal B meson mass j28[, the correct B is chosen in 
75% of all events with multiple candidates. If both possible dipion combinations of the four 
pions fall within the fit region, the dipion with the largest momentum difference between 
its pions is chosen. Since our selection criteria biases the M^ c distributions, we do not use 
this variable in the fit to data. The fraction of mis-reconstructed signal events is found to 
be 10%. 



We also perform a vertex fit and employ the flavor tagging routine described in Ref. [29 
The tagging information is represented by two parameters, the B^ ag flavor q — ±1, and 
the tagging quality r. The parameter r is an event-by-event, MC determined flavor-tagging 
dilution factor that ranges from r = for no flavor discrimination to r = 1 for unambiguous 
flavor assignment. The data is divided into 7 r-bins, denoted as /. 

We study backgrounds coming from b — > uud channels with the same 
final state as signal, namely B° — > a^ir T , afTr T , b^TT T , /o(980)p°, 
/o(980)/o(980), p°7r + 7r~, / (980)7r + 7r~, 7r + 7r~7r + 7r~. The detection efficiencies e of the con- 
sidered four charged pion final states after applying all mentioned selection criteria are listed 
in Table M 

We also determine correction factors to the efficiency that account for the differences 
between data and MC as calculated by independent internal studies at Belle. The correction 
factors r] arising from particle identification criteria are listed in Table IIII1 

Last, we remove an artificial asymmetry arising mainly from momentum ordering in the 
reconstruction by interchanging the two p° candidates for odd event numbers. 



IV. EVENT MODEL 

The branching fraction is extracted from a six- dimensional extended unbinned maximum 
likelihood fit to AE, J^s/b, Mi, M 2 , Hi and H 2 in the I th r-bin and SVD configuration 
s where and Hk represents the invariant dipion mass m^+^- and helicity parameter 
cos6h, respectively, of the p° candidate, identified by index k = 1,2. If linear correla- 
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TABLE III: Correction factors (77) accounting for the difference between data and MC in particle 
identification. The errors are statistical. For B° — > p°p°, we use the Jl value obtained from the fit 
to data. 



r\ (mode) 


SVD1 


SVD2 


B°- 




0.84 ±0.02 


0.85 ± 0.03 


B°- 


/o/o 


0.84 ±0.02 


0.85 ± 0.03 


B°- 


^foP° 


0.84 ±0.03 


0.86 ±0.03 


B ^ 


/ovr+vr" 


0.84 ±0.02 


0.85 ± 0.03 


B ^ 


p°7r + 7T _ 


0.84 ±0.02 


0.85 ±0.03 


B°^7T 




0.84 ±0.02 


0.85 ± 0.03 



tions between fit variables do not cause a noteworthy bias, the probability density func- 
tion (PDF) for each event i, is taken as the product of individual PDFs for each variable 
V(AE\r s/B ,Ml,M^H[,W 2 ) = V{A&) x V{F S/B ) x V{M{) x V(M\) x V{H\) x V{H\)- 
otherwise, correlations between the fit variables are taken into account as described later 
on. We consider 17 components in the event model, where most resonances are described 
by a relativistic Breit-Wigner 

with a mass-dependent width 

r(m) = r p o( -^-5- ) ( )B 1 {p n + n -,p /fl ); (3) 



Pp° 



m p o 



T p o is the nominal width of the p°, p n + n - {i r n 7T + n -) the momentum (invariant mass) of a 
pion pair and p p o (m p o) the momentum (invariant mass) of the nominal p° coming from 
B° p°p° decays. B\ is a Blatt-Weisskopf form factor as described in Ref. (28| . 

The signal model shape is determined from correctly reconstructed signal MC events. 
The PDF for AE for all four pion final states is taken to be a sum of two asymmetric- 
width (bifurcated) Gaussians incorporating calibration factors that correct for the difference 
between data and MC. These factors calibrate the mean and width of the core bifurcated 
Gaussian. They are determined from a large statistics control sample B° — > D~[K + tt~7t~]tt + 
and are used for all four pion final states. Similarly, the parameters describing the signal's 
core Gaussian are made common among all four charged pion final states. The PDF for 
J^s/b, here and throughout this analysis for all components j, is the sum of two bifurcated 
Gaussians in each r-bin, incorporating calibration factors that correct for the difference 
between data and MC for all BB components as their shapes are determined from MC. The 
shape for all other BB components is fixed from the signal model except for the mean and 
width of their core distribution. The distribution in the Mi-M 2 plane is modelled with a two- 
dimensional (2D) efficiency-corrected product of two relativistic Breit-Wigner functions; the 
distribution in the Hi-H 2 plane is described by Eq. [U corrected by a 2D efficiency histogram 
obtained from fully simulated signal MC. 

The mis-reconstructed model shape is determined from incorrectly reconstructed signal 
MC events and its PDF is a histogram for all fit variables, smoothed if applicable. 
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The continuum model shape was studied with data taken below the T(AS) resonance; 
however its entire shape is free in the fit to extract the branching fraction. The PDF for AE 
is taken to be a first-order Chebyshev polynomial. The m n + n - shape is the sum of a second 
order polynomial, a Breit-Wigner for p° and a Breit-Wigner for / (980). Tsib is modelled 
as described before; in addition, we account for a correlation with m^^-by multiplying its 
mean with the factor a 1 = p Q (Mi+M 2 ) + l, where p is a constant in each r-bin. The PDF for 
Hi and H 2 for continuum and for all other backgrounds are symmetrised 2D histograms. We 
use sets of several histograms in bins of other fit variables if a treatment of the corresponding 
correlation is required. 

The charm BB background shape is determined from a large sample of MC containing 
b — > c transitions and is further divided into a neutral and a charged B sample. For both 
samples, AE is correlated with the helicity angles and therefore its PDFs in each sample is 
formed in different bins of Hi, H 2 . We model the neutral component in four areas of Hi, 
H 2 by first-order Chebyshev polynomials where we also add a Gaussian if \Hi^\ > 0.65 or 
\Hi i2 \ > 0.65 and —0.5 < H 2) i < 0.65 (and vice versa), 

V c ^(AE) HuH2 = (fG(AE,fi,a) + (1 - f){ciCi{AE))) HuH2 , (4) 

and the charged component with the sum of Chebyshev polynomials up to the third order, 

V£?T- {^E) Hl , H2 = ^a(AE) + a(H u H 2 )c 3 C 3 (AE). (5) 

i=l,2 

where the cubic term has the factor a obtained from MC and accounting for the correlation, 
a(Hi,H 2 ) = b(\Hi\ + \H 2 \) if \H 12 \ > 0.5 with a constant b = 0.68 ± 0.01. (Otherwise 
a{Hi,H 2 ) = 1). 

The m n + n - shape of the charm B decays is described by a first order Chebyshev 
polynomial where a Gaussian, accounting for / (600), is added to the neutral decays if 
\H k= i >2 \ > 0.65, 

V^f(m n +^-) HuH2 ee / ffl ,H 2 CiCi(m 7r+7r -)(l - f HliH2 )G(m n + 7T -,n,a). (6) 

The widths of both Gaussians of the Fs/b shapes are multiplied with a constant d l in the 
corners of the Hi, H 2 plane, \Hi )2 \ > 0.7 or H X 2 > 0.5 and H 2 X < —0.5. 

The charmless BB background shape is determined from a large sample of MC events 
containing b — >■ u,d,s transitions (again, divided into a neutral and a charged category). 
The AE distribution of neutral B decays is described by a Gaussian plus a third-order 
Chebyshev polynomial, and that of charged B decays by a first-order Chebyshev polynomial 
pcharmiess = ^c^AE), where ci ee p Mi + p M 2 +p 1 (M 1 x M 2 ) with p = -2.44 ± 0.47 and 
pi = 5.46 ± 1.08. Since the neutral charmless Mi, M 2 distribution is correlated with Hi, 
H 2 , we describe its shape by the sums of Gaussian and Chebyshev polynomials in five bins 
of Hi, H 2 , where in most bins also a correlation between the two 7t + tt~ masses is seen. The 
charged m^^- distribution is modeled by a smooth 2D histogram. 

The PDF shapes for the remaining four-pion states are determined from individually 
generated MC samples. For the decay B° — > p°7r + 7r~, we assume a phase space distri- 
bution and account for this assumption in the systematic uncertainty. Since we include 
mis-reconstructed events in the model, a polynomial is added to their AE PDFs. For the 
non-resonant decay, the correlation with the helicity angles is incorporated by taking a dif- 
ferent width of the core Gaussian for the center of the Hi-H 2 plane. The Mi, M 2 shapes are 
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TABLE IV: List of peaking backgrounds, assumed branching fractions and their expected yields 
for SVD1 and SVD2. 





Mode 




B (xKT 6 ) 


SVD1 


SVD2 


B°- 




ttT 


16.5 ±2.5 


65 


299 


B°- 






1.65 ±1.65 


1 


7 


B°- 


■> bf [tt^t^ 




0.17 ± 1.53 


1 


3 



TABLE V: Summary of yields fixed relative to other yields free in the fit. The mis-reconstructed 
yield is fixed relative to the signal yield, and the charmed and charmless B + B~ background yields 
are fixed relative to their respective B°B° background yields. The central values are obtained from 
MC; the errors are statistically. 



Yield 


SVD1 


SVD2 


/V s 
JV Mis 

t> rcharm,s 

B+B~ 
» Tcharmless,s 

iV B+B- 


(0.09 ± 0.0006)A^ D1 
(1.40 ± 0.07)A^ h o a ™' SVD1 
(0.83 ± 0.08)A^f ess ' SVD1 


(0.09 ± 0.0006)A^ D2 
(1.42 ± 0.04)^ h a ™' SVD2 
(0.86 ± 0.04)A^f CSS ' SVD2 



modeled by 2D histograms, except for B° — >■ f f and B° — > fop , where products of two 
Breit-Wigners are used. Since, for B° — > fop , one p° candidate is a p° and the other a f , 
the correlation between M\ and M 2 is taken into account. The mean of the Gaussian of the 
a\Ti J-'s/b shape is multiplied with a = 1 +pi(M\ + M 2 ), to account for the correlation with 
the masses. 

Besides the dominant four-pion contribution B° — > a^ir T , we also fix the other quasi- 
two-body B decay modes in which one of the decay product is a single pion as listed in 
Table IIV| where we assume the unknown branching fraction of B° — > afn T to be 10% of 



B — > a 1 n T 's. A recent measurement by Belle supports this assumption [13[. The other 
two assumed branching fractions are taken from the current world average in Ref. 28 . 
The total likelihood for 116081 B° — > p°p° candidate events in the fit region is 

p-£^|E is // s Nl > s 

C = H J^T II E NjffT&iAE?, F s/B , M u M 2 , H x , H 2 ), (7) 

l,s ,s ' 1=1 j 

which runs over event i, component j, r-bin I and SVD configuration s. Instead of two 
free signal yields Af| ; for each detector configuration, branching fractions for the free four 
pion final states (j) are chosen as single free parameters and incorporated into the fit with 

N; = B(B° f)N' BB ^, (8) 

where and r] s - are given in Tables [TT] and Ifflt respectively. The fraction of events in each 
r-bin I, for component j, is denoted by /j' s . The fractions of all BB components, 
have been calibrated with the B° — > D~[K + 7r~7T~]7T + control sample. In the fit to data, 
we also float the yields N^, N^^ ,s and N B ^^ ess,s and the parameters of the shape 
of the continuum model. The remaining yields are fixed to the values given in Table |V] as 
determined from MC. 
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V. FIT RESULT 



We perform a six dimensional fit to the data. The projection of the fit results onto AE, 
Mi, M 2 , Hi, H 2 and Tsm are shown in Fig. HI The obtained branching fractions, their 
corresponding yield, upper limits at 90% confidence level are given in Table IVIl together 
with fi of B° p°p°. The statistical correlation coefficients between the observables are 
given in Table IVIII 



TABLE VI: Branching fractions with their corresponding yield, upper limits at 90% confidence 
level and the significances (5) for modes with a positive yield from the fit to data. Furthermore, 
the longitudinal polarization fraction jfx of B° — > p°p° is given. The errors are statistical and 
systematic, respectively, statistical errors are given at the la level. 



Branching fraction (xlO 6 ) 


Events 


UL (xlO" 6 ) 


S(a) 


B(B° -> p°p°) = 1.02 ± 0.30 ± 0.22 
f L = 0.2ll^| ±0.11 


166 


< 1.5 


2.9 



B(B° 
B(B° 
B(B° 
B{B° 
B(B° 



IT 7T 7T 7T 

p°ir + ir~) = 



-3.581^3 ± 2.99 



1.70l|;^ ±5.30 



1.34^-12 ± o. 



/ovr+vr ) x B(/ -> vr+vr ) 
/op°) x B(f -> 7T+7T-) = 0.86 ± 0.27 ± 0.15 
/o/o) x B(/„ ->• vr+vr-) 2 = 0.03 



+0.10 
-0.09 



±0.04 



-25 

33 

-27 

149 

-5 



< 11.7 

< 12.2 

< 3.1 

< 0.2 



< 1 
3.0 



TABLE VII: Statistical correlations between the observables. 



B(B° -> X) 


p°p° h 


7T + 7T 7T + 7T 


p TT + 1T 


/ovr+vr 


/o/o 


foP° 


p°p° 


1 0.45 


-0.03 


-0.34 


0.04 


0.02 


-0.02 


h 


1 


0.06 


0.00 


0.03 


0.01 


-0.08 


7r + 7r _ 7r + 7r _ 




1 


-0.35 


-0.23 


-0.04 


0.00 


p°7r + 7r _ 






1 


-0.07 


0.02 


-0.24 


/ovr+vr" 








1 


-0.36 


-0.20 


/o/o 










1 


0.07 



foP° 1 



The resulting B —> p p branching fraction is quite consistent with the previous Belle 
analysis [2l| and it is also in agreement with the value obtained by the BaBar collabora- 
tion [22|]. The fraction of longitudinal polarization in B° — > p°p° decays is somewhat lower 
than previously measured (differing from the BaBar result by 2.1a) and the branching frac- 
tion of B° — > fop decays is significantly higher than indicated by previous measurements. 
We discuss these results in the next section. For the modes where an upper limit was ob- 
tained, we have improved values for B° — > 7t + 7t~7t + tt~ and B° — > f n + 7T~ than the current 
world averages (28| . 

The relative contributions of the components modeled in this analysis are found to be 
0.1% B° -»■ p°p°, 0.1% B° f p°, 92.8% continuum, 6.7% BB background and 0.3% for 
the remaining four pion final states. 
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FIG. 4: (color online) Projections of the fit to the data. The points represent the data and the 
solid curve represents the fit result. The shaded (red) area shows the B° — > p°p° and the (orange) 
long dashed curve the B° — > /op contribution. The (cyan) dash-dotted curve shows all four pion 
final states, the (dark green) dotted curve shows the non-peaking BB contribution and the (bright 
green) dashed curve shows the total non-peaking background, (a) shows the AE projection in the 
Fs/b > 0.5 region, (d) shows the Fs/b projection in the \AE\ < 0.02 GeV region. Here, the 
dashed curve shows the continuum background contribution. Plots (b,c) and (e,f) show the m 1T + 7T - 
and cosOh projections, respectively, in the T S j B > 0.5 and \ AE\ < 0.02 GeV region. 
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VI. VALIDITY CHECKS AND SIGNIFICANCE 



We investigate the stability of the fit result by removing the components B — > 
7r + 7r~7r + 7T~, fo7T + 7r~, and fofo from the fit, as for those the nominal fit gives negative yields 
consistent with zero. The obtained results are given in Table IVIIII where the fit result 
remains stable. 

TABLE VIII: Fit results from fit to data with the branching fractions of B° — > 7r + 7r~7r + 7r~, 
/o7r + 7r~, and /o/o are fixed to zero. 



Branching fraction (x 10 6 ) 


Events 


B(B° -> p°p°) = 1.05 ± 0.31 (stat) 


169 


f L = 0.23 ± 0.20 (stat) 




B(B° -> p°7r + 7r-) = -0.32 ± 3.00 (stat) 


-6 


B(B° fop ) x B(f -»• tt + tt~) = 0.76 ± 0.28 (stat) 


132 



We repeat the fit while also fixing the B° — > p Q ix + ix yield to zero; this fit result is shown 
in Table ITXl Again, the result is in very good agreement with the nominal fit. 

TABLE IX: Fit results from fit to data with the branching fractions of B° — > tt + it~it + it~ , foir + ir~ , 
fofo and p°n + 7r~ are fixed to zero. 

Branching fraction (xlO -6 ) Events 

169 

132 



Since the branching fraction of B° — > fop° is larger than indicated by the previous 
measurement, we investigate the impact on B° — > p°p° when setting B(B° — > fop ) = and 
vice versa. We obtain consistent results, given in Table |Xl in order to visualize each mode 
separately, the nominal fit projections into the p°p° and p°/o windows are shown in Fig. [51 

TABLE X: Fit results from fit to data with the branching fraction of B° — > fop or B° — > p°p° 
fixed to zero. 



Branching fraction (x 10 6 ) 


Events 


B{B° p°p )\ fop o =0 = 1.01 ± 0.31 (stat) 


162 


f L \ fop o=o = 0.26 ±0.21 (stat) 




B(B° /op°)| p o p o =0 = 0.73 ± 0.21 (stat) 


127 



Also, to measure the potential impact of the choice of the starting parameters, we perform 
an additional likelihood scan with random guesses of the initial parameters, where all fits 
converge to the same result. 

Next we refit the data while fixing fi to either or 1. We obtain B(B° — > p°p°)/ i= o = 
(0.79 ± 0.22) x 10~ 6 and B(B° p°p °) /l=1 = (0.22 ± 0.24) x 10" 6 , respectively, which are 
in agreement with our baseline fit result. 



B(B° -> p°p°) = 1.05 ± 0.28 (stat) 
Sl = 0.23 ±0.19 (stat) 

B(B° -> fop ) x B(f -»• iT+iT-) = 0.76 ± 0.25 
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FIG. 5: (color online) In order to better visualize each mode separately, projections of the fit to 
the data into a p p window 0.6 GeV/c 2 < m^- < 0.88 GeV/c 2 (upper row) and a fop window 
0.6 GeV/c 2 < Mi )2 < 0.88 GeV/c 2 and 0.88 GeV/c 2 < M 2 ,i < 1-08 GeV/c 2 (lower row). The 
points represent the data and the solid curve represents the fit result. The shaded (red) area shows 
the B° — > p°p° and the (orange) long dashed curve the B° — > fop contribution. The (cyan) dash- 
dotted curve shows all four pion final states, the (dark green) dotted curve shows the non-peaking 
BB contribution and the (bright green) dashed curve shows the total non-peaking background. 
Plots (a) and (d) show the projection into AE for J~s/B > 0-5 and (b,c) and (e,f) show into cosOh 
for F s/B > 0.5 and \AE\ < 0.02 GeV). 



We evaluate the statistical significance Sq of the result by taking the ratio of the likelihood 
of the nominal fit (C max ) and of the fit with the signal yield fixed to zero (Co); 




The statistical significances of the B° — > p°p° and B° — y f p° yields are found to be 4.6er and 
3.6cr, respectively. In addition, we demonstrate that the minimum is global by performing 
likelihood scans of all measured four pion final states as well as Jl- Likelihood scans of 
B° — > p°p°, fi and B° — > fop , convolved with the corresponding systematic uncertainties, 
are shown in Fig. El and are used to obtain the total significance of these two modes, which 
are 2.9a and 3.0cx respectively, including all systematic uncertainties. 
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VII. SYSTEMATIC UNCERTAINTIES 

Systematic errors from various sources are considered and estimated with independent 
internal studies and cross-checks. These are summarized in Table IXII This includes the 
uncertainty on the number of produced BB events in the data sample. Contributions to 
the uncertainty in the selection efficiency due to particle identification and tracking are 
determined by independent studies at Belle. 

The uncertainty in the p° shape is determined by varying the fixed mass within its world 
average uncertainty {28[. We account for a difference in the fraction of mis-reconstructed 
events between data and MC by varying this parameter by 5% and repeating the fit. Vari- 
ations in the parametric model shape due to limited statistics are accounted for by varying 
each parameter within their errors. Uncertainties in the non-parametric shapes are obtained 
by varying the contents of the histogram bins within ±ler. The systematic uncertainty due 
to fixing the peaking background yields of B° afir T and B° — > bfir T are estimated 
by varying the branching fraction by its world average error and repeating the fit. For 
B° — > afir T , where only an upper limit is known, the variation is taken as half of the upper 
limit. The fit bias was determined from pseudo-experiments by searching for a difference 
between the generated and fitted physics parameters. For the mode B° — > f p° we find an 
non-negligible bias of +16% which is accounted for by subtracting 8% from the fit result and 
assigning an 8% uncertainty. All other biases are found to be small compared to the statisti- 
cal uncertainty and are therefore treated fully as systematic uncertainties. Furthermore, we 
performed an ensemble test where we replaced B° — > p°tt + 7t~ helicity PDF with one where 
the p° is either longitudinally or transversely polarized to generated MC sets according to 
the fit result. Modes in agreement with zero events, except B° — > p°7r + 7r~, were not gener- 
ated, but left free in the fits. The maximal deviation from the nominal model is taken as the 
uncertainty related to the assumption of the B° — > p°7r + ir~ helicity dependency. Finally, 
the uncertainty from neglecting interference between the four pion final states was estimated 
by constructing a 4-body amplitude including detector effects and generating three relative 
interference configurations between two of the four pion final states: maximum constructive 
interference, no interference and maximal destructive interference. We consider the modes 
B° — > p°p°, ai7r, f p° and p°7m where for p°p° we set /x = 0.21 and generate MC sets for 
each combination. For each mode, the largest deviation from the sample with no interference 



17 



TABLE XI: Systematic uncertainties of the branching fraction of B — > p p and /l- 



Category 


-> p°p°) (%) 


SfL 


N(BB) 


1.4 




Tracking 


1.4 





Particle identification 


2.5 





Mis-reconstruction fraction 


1.3 


< 0.001 


p° shape 


0.2 


< 0.001 


Model shape 


5.1 


0.08 


Histogram shape 


5.2 


0.03 


B{B° aivr) 


0.4 


0.03 


B{B° -> 6itt) 


<0.1 


< 0.001 


5(5° -> a 2 vr) 


<0.1 


< 0.001 


Fit bias 


1.9 


0.03 


Interference 


19.2 


0.03 


p°inr helicity 


6.3 


0.05 


Total 


22.0 


0.11 



when fitting all sets with a fitter that assumes incoherent decays, gives the systematic uncer- 
tainty from interference. We found that interference with a\% gives the largest uncertainty 
in all cases. 



TABLE XII: Systematic uncertainties (%) for the peaking background modes. 



Category 


47T 


p°7r + 7T 


/ovr+vr" 


/o/o 


fo P ° 


N(BB) 


1.4 


1.4 


1.4 


1.4 


1.4 


Tracking 


1.4 


1.4 


1.4 


1.4 


1.4 


PID 


2.5 


2.5 


2.5 


2.5 


2.5 


Fit bias 


19 


10 


7 


100 


8.5 


Model shape 


28.5 


218.8 


13.5 


13.8 


6.5 


Histogram shape 


38.1 


127.3 


54.5 


46.1 


5.9 


B(B° -> am) 


10.3 


129.5 


3.1 


4.7 


3.4 


B(B° -> &itt) 


<1 


1.6 


<1 


<1 


<1 


B(B° a 2 Tr) 


<1 


2.7 


<1 


<1 


<1 


Interference 




93.16 






6.43 


p°7T7r helicity 


26.8 


23.3 


17.6 


14.1 


4.5 


Total 


83.5 


312.3 


73.07 


116.9 


15.9 
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VIII. 2 CONSTRAINT 



We use the longitudinally polarized fraction of our result B(B° — > p°p°)\f L= i = (0.21 ± 
0.36) x 10~ 6 to obtain a new constraint on the CKM angle 02 through an isospin analysis (l6| . 
For the remaining sides of the triangle we use the current world averages from the Heavy 
Flavor Averaging Group, Ref. j30[; the longitudinally polarized fraction of B(B° — > p + p~) = 
(2.42 ± 0.31) x 10~ 5 with /+" = 0.978±g;g||, A%~ P = 0.06 ± 0.13, S+p = -0.05 ± 0.17 and 
B(B ± p ± p°) = (2.4 ±0.2) x 10" 5 . The CP asymmetries of B° p°p° are taken from the 
BaBar measurement [l2j; A°£ P = -0.2 ± 0.85, <S£° P = 0.3 ± 0.73. Fig. [7] shows a probability 
scan of 02 that is consistent with the SM in its one a interval of 02 = (91.0 ± 7.2)°. The size 
of the penguin contributions are negligible: A02 = (0.0 ± 5.4)°. 




FIG. 7: 1— CL vs. 02 obtained from an isospin analysis from B — > pp decays. The horizontal line 
shows the 66% CL. 



IX. CONCLUSION 

We have presented a measurement of the branching fraction of B° — > p°p° decays and 
the fraction of longitudinally polarized p mesons in this decay, together with the other four 
pion final states B° p°7r + 7r~, B° f p°, B° f f , B° / 7r + 7r~and B° u+ij-it+it- 
using the final Belle data set of 772 x 10 6 x BB pairs. All results are listed in Table N\\ 
where we have upper limits except for B° — > fop - There, we report the first evidence with 
a significance of 3.0 standard deviations. We used the result of longitudinally polarized p 
mesons in B° — > p°p° decays to further constrain the CKM angle 02 = (91.0 ± 7.2)° with 
an isospin analysis. This measurement is still statistically limited and more insight on the 
interesting and complex structure of four pion final state B decays will be accessible by 
future experiments, e.g, by performing a four-body Dalitz analysis. 
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